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Abstract. - We show that in presence of an applied external field the two-component order 
parameter superconductor falls in two categories of ground states, namely, in the traditional 
Abrikosov ground state or in a new ground state fitted to describe a superconducting layer with 
texture, that is, patched regions separated by a phase difference of 7r. The existence of these two 
kinds of ground states follows from the sole assumption that the total supercurrent is the sum of 
the two individual supercurrents and is independent of any consideration about the free energy 
expansion. Uniquely defined relations between the current density and the superfluid density hold 
for these two ground states, which also determine the magnetization in terms of average values 
of the order parameters. Because these ground state conditions are also Bogomolny equations 
we construct the free energy for the two-component superconductor which admits the Bogomolny 
solution at a special coupling value. 
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Introduction. — Many thermodynamic and trans- 
port phenomena in superconductivity can be explained on 
the basis of a macroscopic order parameter, as shown by 
the phenomenological approach of Ginzburg and Landau 
(GL), proposed many year before Leon Cooper found that 
a weak attraction binds electrons near to the Fermi sur- 
face to form pairs. From the GL theory A. A. Abrikosov 
predicted the existence of two types of superconductors 
according to the dimensionless coupling k, thus proving 
that this theory is well suited to describe alloys, which are 
k > 1/V2 superconductors, able to sustain a stable vortex 
state in presence of an external applied field [T]. Central 
to the present analysis is the observation that at the core 
of Abrikosov's treatment lives a condition which is totally 
independent of the Ginzburg-Landau free energy. This 
condition is the requirement that the macroscopic order 
parameter fulfils a ground state equation and from it the 
supercurrent and the superfluid densities are shown to be 
directly related to each other. It also follows from this 
condition that the magnetization is proportional to the 
spatial average of the superfluid density, a mean field re- 
sult that remains valid in the presence of thermal fluctua- 



tions [2j . Therefore the ground state condition for the one- 
component order parameter (1COP) superconductor lives 
at a more fundamental level than the Ginzburg-Landau 
free energy expansion. To attain full predictive power from 
the magnetization expression the superfluid density must 
be determined in terms of fundamental parameters, and, 
at this point, the free energy considerations becomes valu- 
able, as shown by Abrikosov, who obtained from the GL 
theory the magnetization as a function of the temperature 
T and of the external applied magnetic field, H, taken to 
be near to the upper critical field H C 2 [3]. 

In this letter we show that for the two-component or- 
der parameter (2COP) superconductor the requirement 
that the macroscopic order parameter fulfils a ground 
state equation leads to just two possibilities, namely, the 
Abrikosov ground state, applied individually to each one 
of the components, and a new ground state that describes 
a textured superconducting layer because of the intrinsic 
phase difference of ir between its distinct regions. There- 
fore we find here that for the 2COP superconductor the 
ground state condition also exists and lives in a more fun- 
damental level than the free energy expansion. In case 
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the local field is equal to the applied one, the Abrikosov 
ground state locks the 1COP into the lowest Landau level, 
whereas this new 2COP ground state takes contributions 
from all Landau levels. 

In 1976 E. Bogomolny [3HZ] found, while working in 
string theory, extra properties for the GL theory valid for 
K = l/y/2: first order equations, instead of the second or- 
der variational equations (Ampere's law and the GL equa- 
tion) solve the GL theory. The Bogomolny's first order 
equations share a common feature with the Abrikosov's 
treatment of the GL theory: one of these equations is ex- 
actly the Abrikosov ground state condition. 

We also study here a 2COP free energy because the 
new ground state condition is one of its Bogomolny equa- 
tions. This means that at a particular coupling, defined 
to be k = l/y/2, like for the 1COP case, the minimum of 
this 2COP free energy is reached by first order equations 
instead of second order ones, which are then defined as 
2COP Bogomolny equations. Therefore this 2COP free 
energy theory is a unique generalization of the 1COP GL 
theory. Previously proposed 2COP GL theories [8HT2"] do 
not contain this feature, namely, none of them is mini- 
mized by the present Bogolmony equations associated to 
this new ground state. We stress that the study of this 
2COP free energy is just complementary to our major 
claim of a ground state condition independent of a free 
energy proposal. 

Nowadays a broad range of experimental data is point- 
ing towards common features for the superconductors 
whose pairing originates in two-dimensional layers, such 
as the existence of two fundamental energy scales [T31H4] 
and the competition of the superconducting state with a 
magnetic collective state [T5]. This leads to akin proper- 
ties for the layered superconductors albeit their distinct 
atomic structure. The most studied materials are the 
cuprates whose two-dimensional layers are made of copper 
and oxygen atoms [T6HT8] . Recently the pnictides brought 
excitement into the field because superconductivity origi- 
nates in two-dimensional layers made of iron and arsenic 
atoms instead [H3H2T] . The origin of superconductivity, 
the microscopic nature of the pairing mechanism [2"Tll2"8] . 
and the intrinsic magnetism are under current investiga- 
tion [HQSHS]- We propose here that the new 2COP 
ground state will help to understand the macroscopic 
properties of the cuprates and of the pnictides because 
of its intrinsic layered structure. 

The new ground state is described by a 2COP, ^ = 
(| ip\ | exp (iOi), \ipn \ exp (162)), whose response to the pres- 
ence of an applied field H = Hx 3 is a local magnetization 
proportional to ^ . In components this magnetization 
is given by, 

M\ = -(ft?/mc)(cos^i||^|), (1) 
M 2 = (ng/mc)(sinfl|Vi||^|>, and, (2) 
M 3 = -(%/2mc)(|^| 2 - \H 2 ), (3) 

where 8 = 6\ — 82 is the phase difference between the two 



order parameters and (• • •) means spatial average. There- 
fore there is a local transverse magnetization to the applied 
field, but the total average magnetization should vanish, 
Mi = M2 = 0, to avoid a spontaneous torque in the sys- 
tem. This is a novel feature, not present in the Abrikosov 
ground state, whose magnetization is fully oriented along 
the applied field. The key ingredient to make the aver- 
age in plane magnetization to vanish is to introduce tex- 
ture, which means that the local transverse magnetization 
flips sign by the phase shift 8 — >• 8 + ir. Phase texture in 
two-gap superconductors has been discussed in the liter- 
ature [52]. It introduces interesting features such as the 
onset of solitons [2lfll2"4"] and an intercomponent Josephson 
interaction |25j that may directly relate to an abnormal 
AC loss peak [26] . 

The 1COP GL theory has been applied to the cuprates 
cither by including a mass anisotropy or by the Lawrence- 
Doniach model, which treats superconducting layers cou- 
pled by the Josephson effect pQ. The present 2COP free 
energy allows for such generalizations that will be consid- 
ered elsewhere. 

One-component order parameter (1COP) 
ground state. — On the basis that the superconduct- 
ing state can be described by a macroscopic wavefunction, 
ip, the supercurrent density is J = (q/2m) 

D = (h/i)W - (q/c)A. Without invoking the GL theory, 
Ampere's law, V x h = (4tt/c)J, is exactly solved for the 
local field, h — V x A, in terms of -0 by assuming the 
ground state condition, 

D+^j = 0, (4) 

where D± = Di±iD2- This ground state condition applies 
for a bulk superconductor with continuous axial symmetry 
along the applied field direction (H = Hx^), such that 
fields only dependent on the coordinates orthogonal to it, 
(xi,X2). From the real and imaginary parts of Eq.Q, 
we find that the supercurrent and superfluid densities are 
related by, 

J = "^V x (M 2 x 3 ) (5) 
Integration of Ampere's law gives that the local field is, 

^ 3 = i?-— |Vf, and 4irM = — — (H 2 ). (6) 

mc mc 

is the magnetization. For the 1COP case the spatial av- 
erage value is over the orthogonal coordinates, such as for 
the magnetic induction, B = H + 4irM, which becomes 
B = (h 3 )x 3 = (J d 2 xh 3 / S)x3, since along the third di- 
rection and under translational invariancc, it suffices to 
consider integration over the unit cell area S. For a con- 
stant field (Ai = —Hx2, A 2 = 0), Eq.dH is just the lowest 
Landau level condition whose solution is 

^^c^-^+w) 2 . (7) 

k 
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The set of wavenumbers k and the constants Ck are de- 
termined by imposing periodic conditions to the order pa- 
rameter and fixing the number of vortices within the unit 
cell area. One obtains that B = N&o/S, for N vortices 
within the cell. Notice that the order parameter of Eq.© 
together with the local magnetic field of Eq.© stem just 
from the ground state condition. Nevertheless they de- 
scribe the vortex state without invoking the free energy, 
that only enters to select the vortex state of lowest free en- 
ergy. As shown by E.H. Brandt pQ, the order parameter 
of Eq. ([7]) and the local field of Eq. © provide an excellent 
description of the full GL free energy solution for fields 
in the range Q.5H c2 < H < H C 2- Therefore the ground 
state condition of Eq.Q does describe the GL theory for 
an applied field not necessarily near to H C 2- 

Two-components order parameter (2COP) 
ground state. — Along the crystal's major axis, where 
the mass tensor is diagonal, and an appropriate combina- 
tion of the two order parameters is taken, the supercurrent 
is expressed by J = J\ + J 2 , Jj = (q/2m) ypjD^j + c.cj , 
j = 1, 2. The Abrikosov ground state assumes trans- 
lational invariance along the applied field such that the 
condition of Eq.© applies to both components, and the 
local field becomes h 3 = H — (hq/ mc) ( | tpi | 2 + | tp2 | 2 ) ■ This 
new ground state satisfies the condition, 



a ■ = 0, where * 
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where a are the Pauli matrices, the generators of 
the SU(2) group: ai=( ® J, a 2 =( °. M, and 



CT3 



These matrices live on a vector repre- 
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scntation of the 0(3) group, and so this ground state 
contains an orientation and a direction in space, asso- 
ciated to the physical superconducting layer linked to 
the crystal's major axes. Take the two components of 
Eq.©, D 3 V»i + D_i> 2 = and £>+Vi - D 3 ip 2 = 
to separate the two order parameters in second order 
equations: [(D? + D%) + D\ - (qh/c)h 3 ]tpi = and 



[(Dl+D$)+D% + (q%/c)ha]il)2 = 0. We show that this new 
2COP ground state is truly three-dimensional, although it 
describes fields that evanesce from a layer, opposite to the 
Abrikosov ground state, which is two-dimensional and de- 
scribes a translational invariant bulk state along the third 
direction. We obtain the exponential decay away from the 
layer in presence of H such that the local field is domi- 
nated by the applied field, h 3 ps H > 0. The kinetic energy 
within the orthogonal plane is positive, (D\ + D 2 )ipi > 0, 
i = 1,2, and so, we must have that D^tp 2 < in order to 
obtain a non-trivial solution for the second equation. We 
take the gauge A 3 = and that the layer is at x 3 = 0, 
such that tpi — cxp (— \q 3 \x 3 )^i(xi, x 2 ), i = 1,2, to obtain 
that, 



(D( + Di) in 



(9) 



(Dl + D\)4>2 = [(hq 3 ) 2 -jh 3 ]4j2, (10) 

For the case h 3 = H the system admits the general solu- 
tion, 

h W(n-i)(«i,«a) J' 1 11J 



he ' 1 ^ hq 3 (n) 



D+ijj {n) , (12) 



where (D\+ Dl) 1 ^ ( n ) = (hqH / c)(2n+l)ip^ and V'(-i) = 
by definition. The first term retrieves the lowest Landau 
level condition for the 1COP (q 3 = 0, ip2 = and D+V-'i = 
0) and so the above 2COP contains both two and three 
dimensional contributions. We stress that each pair (tp( n ), 
V>(n-i)) is bound by the same coefficients Ck and wave 
numbers k, but not necessarily two different pairs. 

The major result of this letter is that under the ground 
state condition of Eq. ((5J , the supercurrent becomes 



2m v ; 



(13) 



Integration of Ampere's law gives that, 

h = H- ^*t^ ail d A-kM = -^(^u^),(U) 

vac ' mc \ / ' V / 



for the local magnetic field and for the magnetization, re- 
spectively. The latter is a three-dimensional average value 
defined as B = (h) = J d 3 xh/V for the magnetic induc- 
tion, where V is the volume of the unit cell. 

In conclusion several ground states with elaborate mag- 
netic patterns can be constructed from our explicit 2COP 
solution ^E", given by Eqs.([TT]) and (TT21 . As previously dis- 
cussed these ground states are textured which means that 
they have neighbor regions phase separated by 7r in or- 
der to render that /*t CTl ^^ = /^rt^^rN = o, yielding 

a vanishing magnetization along the layer. This in plane 
magnetization must vanish inside the bulk of the super- 
conductor, and there are several ways to do it, either by 
considering two distinct layers with opposite in plane mag- 
netization, stripes [27] or checkerboard patterns [29] inside 
a single layer. 

Virial relation. — Wc provide another derivation of 
Eqs.([6]) and (|T4]) from the ground state conditions, Eqs.Q 
and ([5]). This derivation brings insight into the kinetic en- 
ergy and is based on the scalar virial relation [3"0H3"2"] in- 
stead of the supercurrent density. The scalar virial relation 
holds for both the 1COP and the 2COP and states that 
the applied field times the magnetic induction is propor- 
tional to the average kinetic energy plus twice the average 
field energy: 



HB 

An 



2m 



8tt 



(15) 
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The kinetic energy is expressed along the crystal's ma- 
jor axis, taken appropriate combinations of the two com- 
ponents, that bring it to the above form. We briefly 
review the 1COP case, considered by U. Klein and B. 
Pottinger |33j sometime ago. The virial relation is 
HB/Att = (IL'iV'I 2 + \D 2 ijj\ 2 )/2m + (h 2 )/4ir, and from 
it, we obtain Eq.©. To show this, express the ki- 
netic energy in terms of the ground state condition of 
Eq.Q: eJ=i \ D ^\ 2 = E?=iVj (irDjfp) + \D+T/A 2 + 
(hq/c)h 3 \ip\ 2 . The total derivative term does not con- 
tribute to the average kinetic energy because of peri- 
odic boundary conditions. One obtains that H(h 3 )/Air = 
(hq/2mc)(h 3 \ip\ 2 ) + (ti 1 ) /Air. Then one easily verifies that 
the solution for this relation is indeed given by Eq.©. 
Similarly the average kinetic energy of Eg . ([15)) can be ex- 
pressed through the ground state condition, using that, 

(\D^\ 2 ^ = (\ff- D^\ 2 ) + (hq/c)(h • tftfftfV (16) 

The virial relation of Ea. (|T5j) leads to the local field of 
Eq . (fT4")l , by introducing Eq.© into the above expression. 

The kinetic energy density is intrinsically positive, and 
so, helpful to understand why in presence of H the 2COP 
has the top component in a Landau level higher than the 
bottom one, according to Ea.(fTTj). Take a constant ap- 
plied field along the third direction in Eq.©, use Eqs. (fl~4l 
and (|16p , and neglect the fourth order term under the as- 
sumption of weak order parameter because of the proxim- 
ity to the normal state. It follows that (\D^/\ 2 ) /(Hq/c) « 
#(|?/>i| 2 - \tp 2 1 2 ), indicating that (\ipi\ 2 ) > (|V>2 1 2 )- An in- 
teresting property is that {^a^ty) = holds for each pair 
* given by (?/>(„), D + ip {n) /ihq 3 (n)). 

Ginzburg-Landau Theory. — The 2COP GL theory 
with Bogomolny solution has three temperature indepen- 
dent parameters, oti(T), i = 1,2,3, that form a vector 
and consequently introduce intrinsically preferred direc- 
tions for the superconductor. Nevertheless it has a single 
dimensionless coupling that for k = l/\/2 admits a Bogo- 
molny solution, as shown below. 

We briefly review the derivation of the Bogomolny so- 
lution for the 1COP GL theory, whose free energy is 
F(B, T) = (\DiP\ 2 /2m ~ a{TM 2 + f3\i>\ 4 /2 + h 2 /8n), 
where a(T) = (T c — T) near to the transition and 
P > 0, rendering possible to have a non-trivial solu- 
tion ip in case that u(T) > 0. Expressing quantities 
in dimensionless units based on the coherence length, 
£(T) 2 = h 2 /2m\a(T)\, and the thermodynamic magnetic 
field, H c (T) 2 /4ir = a(T) 2 /^, one finds that the GL 
theory only depends on a single dimensionless temper- 
ature independent coupling k = H c2 (T) / {^/2H C (T)) = 
(mc/hq)yj fi/{2ir), the upper critical field being H C 2(T) = 
$o/2ir£(T) 2 . Hereafter we use the dimensionless quan- 
tities: x -> xji (V -> £V), A -> (2ir£/<Z> )A (D -> 
(h/QD, D -» V/i — A, h -> h/H c2 (H -> H/H c2 ), 
F -> F/(H 2 /4tt) (G -> G/(H 2 /4ir)). For instance, the 
thermodynamic relation H = AndF/dB becomes in re- 



duced units H = (l/2n 2 )dF/dB. Consider the special 
case of symmetry along the third direction defined by H 
and decompose the kinetic energy using the key condi- 
tion to obtain that the free energy is F = (\D + ip\ 2 + 
MV'I 2 — IV'I 2 + IV'| 4 /2 + K 2 h 2 ). A straightforward rear- 
rangement of the terms leads to F = —1/2 + B + Fq, where 
^0 = (\D + i;\ 2 ) + ((h 3 + H 2 - l) 2 )/2 + (^-l/ml)- For 
k = 1/ \[2 Fq is positive and hold the first order differen- 
tial equations, D + ip = and h 3 + \ip\ 2 — 1 = 0, at the 
minimum since Fq = 0. Consequently, at this particular 
coupling, vortices do not interact as the free energy is ad- 
ditive in their number, N, and F = —1/2 + B, where the 
magnetic induction is B = NQq/S. 

The 2COP GL theory which has a Bogomolny symmetry 

is, 




The out of plane coupling, a 3 , is associated to T c , 
whereas the in plane ones, ct\ and a 2 , are related to 9 
transitions. The simplest of all situations, a\ — a 2 = 0, 
brings some insight into the problem since the remaining 
quadratic coupling becomes — a3(T) (|i/>i| 2 — |i/> 2 | 2 ) . In 
this case ip± admits a non-trivial solution below T c 
(a 3 (T) > 0), whereas i\) 2 does not, and must vanish. 
However above T c (a 3 (T) < 0) a non-trivial solution 
for ip 2 must be prevented and this is done by assuming 
that a 3 vanishes above T c . We get further insight into 
the properties of this theory by doing a SU(2) rotation 
= C/(T)tf, such that a(T) ■ U(T)aU(T)^ = \a(T)\a 3 . 
Under this transformation the free energy becomes 
F = (\m R \ 2 /2m- \d{T)\^ R a^ R + /3(^ R ) 2 /2+ 
+h 2 /8ir). According to the previous argument we 
conclude that no stable second component is possible, 
namely, ipR 2 = turning this theory essentially into a 
1COP GL theory for ip Rl . Notice however an important 
difference, does not vanish at T c because of the 

presence of a\ and a 2 , assuming that these parameter 
only vanish above T c . The previously defined quantities 
HJT), £(T), and H c2 (T) depend on \a(T)\ instead and 
this leads to the same dimensionless n of the 1COP 
case. For instance the one-dimensional surface energy 
barrier for the 2COP theory of Ea. (fl7]) is found to be 
a independent and only k dependent, exactly like for 
the 1COP case, because the above rotation is possible. 
However for the three-dimensional problem we demand 
from the theory that the ground state condition of Eq.® 
applies to ^ and not to ^ R . This will halt the rotation 
preventing the theory to become a one component free 
energy, but keeping its single scale properties. Reduced 
units are defined as before, and the free energy becomes, 

F = ^|l3*| 2 - yfy • *to^ + 1 (ijt^) 2 + K 2p\ (18) 

For zero magnetic field and constant order parameter the 
minimum of the 2COP GL potential gives the values for 
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the two physically relevant order parameters. Using that 
(* t *) 2 = (\fdo i \f f ) 2 , the GL potential becomes, F P = 
— 1/2 + (^(T^ — a/\a\) 2 /2, whose minimum is reached at 
|Vi| 2 = (H + a 3 )/2|a|), \H 2 = (l«l " a 3 )/2\S\) and 
tan 6 = —a.<2,l<y.\. Notice that it is always true that > 
| ip2 1, and that above T c they become equal, assuming that 
ol\ and oli exist above the transition. This restores a U(l) 
invariance to the theory. A change of sign in a\{T) or 
&2(T) changes the phase 9 but not the minimum of \tpi\ 
and |t/>2 I- The transformation oi\ — > — ol\ and a2 — > —0:2 
is a symmetry of the minimum, equivalent to a — > 9 + ir 
rotation. This theory has a Bogomolny solution for k = 
1/V2, given by a ■ = and h + ^a^ - a/\a\ = 0, as 
the free energy can be expressed by, 

F = (\<f-D^\ 2 + [h+ (^a^-d{T)/\d{T)\)] 2 /2' / 

+ (k 2 - 1/2) (h 2 ) + (a(T)/\a(T)\) ■ B/2 - 1/2. (19) 

Thus the 2COP free energy studied here is distinct from all 
other free energies studied so far [5HT2"] because the ground 
state condition holds in two distinct applied field H and 
k regimes, showing that it is a unique generalization of 
the 1COP GL theory: k > 1/^2 and 0.5H c2 < H < H c2 
(Abrikosov), and also n = l/v2 throughout the whole H 
regime (Bogomolny). 

In conclusion, wc have found here a textured ground 
state for the two-component order parameter whose in- 
trinsic transverse magnetic moment averages to zero along 
the superconducting layer because of 7r phase difference 
between distinct regions. This ground state solves one of 
the Bogomolny equations of a free energy obtained here. 
The ground state condition lives in a more fundamental 
level than that of the free energy expansion. 



We thank Ernst Helmut Brand and Said Salem Sugui Jr. 
for helpful discussions. Part of this work was supported by 
the bilateral project between Flanders and Brazil. M.M. 
Doria also thanks CNPq, Capes and FAPERJ (Brazil). 

REFERENCES 

[1] Ernst Helmut Brandt, Rep. Prog. Phys., 58 (1995) 
1465. 

[2] Baruch Rosenstein and Dingping Li, Rev. Mod. Phys., 
82 (2010) 109. 

[3] P. G. de Gennes, Superconductivity in Metals and Alloys 

(Persus Book,Massachusetts U.S.A.) 1989. 
[4] E. B. Bogomolny, Sov. J. Nucl. Phys., 24 (1976) 449. 
[5] Laurence Jacobs and Claudio Rebbi, Phys. Rev. B, 19 

(1979) 4486. 

[6] I. Luk'yanchuk, Phys. Rev. B, 63 (2001) 174504. 

[7] F. Mohamed, M. Troyer, G. Blatter and I. 
Luk'yanchuk, Phys. Rev. B, 65 (2002) 224504. 

[8] M. P. Das, Hong-Xing He and T. C. Choy, Interna- 
tional Journal of Modern Physics B, 2 (1988) 1513. 



[9] 
[10 

[11 

[12 

[13 
[14 

[15 
[16 

[ir 

[18 

[19 
[20 



[21 
[22 

[23 
[24 

[25 

[26 

[-" 

[28 
[29 

[30 

[31 

[32 

[33 



Manfred Sigrist and Kazuo Ueda, Rev. Mod. Phys., 63 
(1991) 239. 

A. J. Berlinsky, A. L. Fetter, M. Franz, C. Kallin, 
and P. I. Soininen, Phys. Rev. Lett, 75 (1995) 2200. 

Joynt Robert and Louis Taillefer, Rev. Mod. Phys., 
74 (2002) 235. 

Egor Babaev and Martin Speight, Phys. Rev. B, 72 

(2005) 180502. 

S. Hufner, M. A. Hossain, A. Damascelli, and G. A. 
Sawatzky, Rep. Prog. Phys., 71 (2008) 062501. 

A. Bussmann-Holder, A. Simon, H. Keller and A. 
R. Bishop, Journal of Superconductivity and Novel Mag- 
netism, 23 (2010) 365. 

Eun Gook Moon and Subir Sachdev, Phys. Rev. B, 80 
(2009) 035117. 

Y. Kohsaka, C. Taylor, K. Fujita, and A. Schmidt, 
C. Lupien, T. Hanaguri, M. Azuma, M. Takano, H. 
Eisaki, H. Takagi, S. Uchida, and J. C. Davis, Science, 
5817 (2007) 315. 

S. D. Wilson, P. Dai, S. Li, S. Chi, H.-J. Kang, and 
J. W. Lynn, Nature, 442 (2006) 59. 

D. Haug, V. Hinkov, A. Suchaneck, D. S. Inosov, 
N. B. Christensen, Ch. Niedermayer, P. Bourges, 
Y. Sidis, J. T. Park, A. Ivanov, C. T. Lin, J. Mesot, 
and B. Keimer, , Phys. Rev. Lett, 103 (2009) 017001. 

T.-M. Chuang, M. P. Allan, Jinho Lee, Yang Xie, 
Ni Ni, S. L. Bud'ko, G. S. Boebinger, P. C. Canfield, 
and J. C. Davis, , Science, 327 (2010) 181. 

A. J. Drew, F. L. Pratt, T. Lancaster, S. 
J.Blundell, P. J. Baker, R. H. Liu, G. Wu, X. H. 
Chen, I. Watanabe, V. K. Malik, A. Dubroka, K. W. 
Kim, M. Rossle, and C. Bernhard, , Phys. Rev. Lett., 
101 (2008) 097010. 

T. Yildirim, Phys. Rev. Lett, 101 (2008) 057010. 

A. Gurevich and V.M. Vinokur, Phys. Rev. Lett., 97 

(2006) 137003. 

Y. TANAKA, Phys. Rev. Lett, 88 (2001) 017002. 
A. Gurevich and V.M. Vinokur, Phys. Rev. Lett., 90 
(2003) 047004. 

Helmut R. Brand, Mauro M. Doria, and Harald 
Pleiner, Phys. Rev. Lett, 60 (1988) 2810. 

Y. Tanaka, A. Crisan, D. D. Shivagan, A. Iyo, K. 
Tokiwa, and T. Watanabe, Jpn. J. Appl. Phys., 46 

(2007) 134. 

Erez Berg, Eduardo Fradkin, Steven A. Kivelson, 
and John M. Tranquada, New Journal of Physics, 11 
(2009) 115004. 

Michael R. Norman, Science, 325 (2009) 1080. 

E. Bascones and B. Valenzuela, Phys. Rev. B, 77 

(2008) 024527. 

M. M. Doria, J. E. Gubernatis, and D. Rainer, Phys. 
Rev. B, 39 (1989) 9573. 

Mauro M. Doria and Sarah C. B.de Andrade, Phys. 
Rev. B, 53 (1996) 3440. 

Mauro M. Doria, Ernst Helmut Brandt, and F. M. 
Peeters, Phys. Rev. B, 78 (2008) 054407. 

U. Klein and B. Pottinger, Phys. Rev. B, 44 (1991) 
7704. 



p-5 



